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Abstract: The electroweak symmetry breaking transition may supply the appropriate 
out-of-equilibrium conditions for baryogenesis if it is triggered sufficiently fast. This can 
happen at the end of low-scale inflation, prompting baryogenesis to occur during tachyonic 
preheating of the Universe, when the potential energy of the inflaton is transfered into 
Standard Model particles. With the proper amount of CP-violation present, the observed 
baryon number asymmetry can be reproduced. Within this framework of Cold Electroweak 
Baryogenesis, we study the dependence of the generated baryon asymmetry on the speed 
of the quenching transition. We find that there is a separation between "fast" and "slow" 
quenches, which can be used to put bounds on the allowed Higgs-inflaton coupling. We 
also clarify the strong Higgs mass dependence of the asymmetry reported in a companion 
paper M. 
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1. Introduction 

The possibility of generating the observed baryon asymmetry at the electroweak phase 
transition dates back just over twenty years [Q]. The Standard Model provides baryon 
number violating processes as well as CP-violation and potentially departure from thermal 
equilibrium. The task at hand is to discover a scenario in which these ingredients come 
together and generate the correct amplitude of the baryon asymmetry. 

Cold Electroweak Baryogenesis was proposed as an alternative to Electroweak Baryo- 
genesis at a bubble wall in a first order phase transition [||, [5j . In the Minimal Standard 
Model such a transition is ruled out by bounds on the Higgs mass || 0, g, and Cold 
Electroweak Baryogenesis relies on an out-of-equilibrium symmetry breaking transition 
triggered by the evolution of an inflaton field. 

Inflation is supported by measurements of the CMB Q and although the precise origin 
of the accelerated expansion has not been confirmed, models based on a scalar field rolling 
in a suitable potential can successfully reproduce observational signals. The energy scale of 
inflation is constrained by observations, but is ususally guided by expectations of physics 
beyond the Standard Model. In Cold Electroweak Baryogenesis, we apply a minimal exten- 
sion to the Standard Model, a single inflaton, and assume that the energy scale of inflation 
is the Electroweak scale. Although this does not discard the possibility of further physics 
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at higher energies, it allows for a self-contained description of baryogenesis. This because 
a period of electroweak-scale inflation will reduce the impact from whatever came before. 

In section 2 we will introduce the relevant details of Cold Electroweak Baryogene- 
sis. Tachyonic preheating with varying quench times is studied in section 3, with no 
CP-violation, where we observe the behaviour of CP-even observables. These results are 
then used in section 4, where we add CP-violation and compare the full simulation results 
to linear treatments based on the CP-even observables, as in jjj. We interpret the quench- 
time and mass dependence in terms of the winding around zeros of the Higgs field. Our 
conclusions are in section 5. 



2. Cold Electroweak Baryogenesis 

A number of conditions are to be met for Cold Electroweak Baryogenesis to be successful: 

• A period of inflation has to take place that ends at sufficiently low energy, such that 
the reheating temperature is well below the electroweak scale. Then baryon number 
changing processes (sphaleron transitions) cannot wash out a previously generated 
asymmetry. Such low scale inflation can be made to agree with observational con- 
straints, with some allowance for parameter tuning |5], [TT|| , 

• Baryon number violation has to be effective while the system is out of equilibrium. 
The finite-temperature transition is a crossover rather than a phase transition in the 
Minimal Standard Model, f], 0, P, and hence the system stays close to equilibrium. 
The phase transition can however be first-order in supersymmetric models with an 
enlarged Higgs sector (see for instance |l2| ). 



We pursue here a different route. Out-of-equilibrium conditions can be introduced by 
triggering electroweak symmetry breaking through a coupling of the Higgs field to a 
rolling inflaton at or after the end of inflation. This results in a tachyonic instabiliy, 
which in turn can result in change of baryon number || f|, ||, 13 , [b|, [j] . 



Although the gauge fields come out of equilibrium, an asymmetry is only created 
in the presence of CP-violation, the strength of which will determine the magni- 
tude of the baryon asymmetry. The quark sector of the Minimal Standard Model 
seems unlikely to provide the required CP-violation Jl5|, [l6|, [l?], |l8|, 19 1. The lepton 



sector has not been studied in similar detail. Here, as in [|l4|, |T|, we use only a sim- 
ple CP-violating interaction composed of Higgs and gauge fields and consider it as 
representing the generic case. 

A simple implementation of the scenario is obtained by adding effective CP violation 
to the SU(2)-Higgs sector of the SM, with action (we use the metric ( — h +- 1-)) 



S = - J d 3 x dt ^ Tr F^F^ + (D^D^ + e + + A 



2 + k^^TiF^F^ 
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The effective mass /x e g is supposed to depend on time, e.g. due to the coupling to an inffaton 
held a. 



VcS = V ~ K<t>o , (2.1) 

with a potential to appropriately generate low-scale inflation. In 'inverted hybrid inflation' 
models g ID), 

a is very small during inflation as it rolls away from the origin. After 
inflation has ended it becomes substantial such that /U 2 g flips sign when a 2 = ^ 2 /A CT ^. 
Much later a and </> settle near their vacuum values and the eigenmodes and eigenvalues of 
the inflaton-Higgs mixing mass matrix determine the masses and other properties of the 



ensuing spinless particles [11] 



The speed of the transition may be characterized by 

(2.2) 



s/2 d\x 2 



off 



nitf dt 



where we used the Higgs mass to set the scale. In [11] we assumed that viable baryogenesis 
would require a sufficiently fast quench \u\ > 0.15. For very fast quenches, the system is 
expected to be very much out-of-equilibrium, and in the limit of infinitely slow quenches, 
the system should stay in equilibrium throughout the symmetry breaking transition. One 
issue is to identify a transition between "fast" and "slow" transitions. This will allow us 
to constrain A^ and/or the shape of the inffaton potential around the time of symmetry 
breaking. Note that in the particular implementation of [O], inflation ends long before 
electroweak symmetry breaking, so that there is some freedom to tune the inflaton speed 
to accommodate a fast quench. 

Setting aside the inflaton and the shape of its potential, we model the effective mass 
parameter by a linear form 

/i 2 ff (t) = /i 2 (l - , < t < tg, (2.3) 

/&(*)= "/A t>t Q . (2.4) 

with tQ introduced as the quench time. In this model the Higgs v.e.v. and mass are given 
by the usual formulas v 2 = fJ> 2 /\ and m 2 H = 2/i 2 , and e = /^ 4 /4A is chosen to set the energy 
density of the vacuum to zero. Furthermore 

1 V2 , v 

u = = — . (2.5) 

fxtg m H t Q 

For \u\ = 0.15, we find rnutQ ~ 9. 

The baryon asymmetry now depends on three parameters 

• The quench time tQ. 



The strength of CP-violation, encoded in the coefficient of the CP-violating term, 

^ss S C p 

3 5cp 



which we write in terms of a dimensionless 5 cp , 



(2.6) 
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In [H| we found that the dependence on S cp is linear at least for 5 cp < 1. For example, 
at ran = 2my/ this leads to 

— = -(0.32 ± 0.04) x 10" 4 8 cp . (2.7) 

n 7 

• The mass of the Higgs boson, written as 

fm H y = 8A 

In P], the mass dependence was seen to be quite complicated in the case of an 
instantaneous quench. We will confirm here, that the mass dependence is also present 
with non-instantaneous quenches; even the sign of the baryon asymmetry depends 
on tuh- 

Shortly after /i 2 ff has turned negative, the occupation numbers of the fields grow (faster 
than) exponentially with time. Initially, the non- linear terms in the equations of motion 
are not very large and it is reasonable to estimate (<jr(t>) using a Gaussian approximation. 
The equation 

4>-V 2 (f> + n 2 eS <t> = 0, (2.9) 



where 



fil s = -M 3 (t-t c ), t c = t Q /2, M = (2.10) 



can be solved analytically, assuming the initial state at t = to be the free-field vacuum 



(see e.g. [21], for the instantaneous quench see p2| ). Following the steps taken in [21], the 
(Fourier transform of the) equal-time two-point function is found to be given by 

(Mc) = lAWI 2 . (2-H) 
f k (t) = C lk Bi(r - k 2 /M 2 ) + C 2k Ai(r - k 2 /M 2 ), t = M(t — t c ), (2.12) 

C lk = ^=\Ai {—Lo\jM 2 ) + iui k Ai{-u 2 k /M 2 )}, u k = V/U 2 + A: 2 , (2.13) 

V 2uJ k 

C 2k = -^[Bi'i-ul/M 2 ) + iLo k Bi(-oo 2 /M 2 )}. (2.14) 

Since the Airy functions behave for large z like Ai(z) cx e^ -2 / 3 ^ 3 , Bi(z) oc e^ 2 / 3 ^ 3 , it 
can be shown that occupation numbers of the unstable modes with k < M^fr grow very 
rapidly once r > 0, and a classical approximation can be made for r > 2 ^T|. I n fact, since 
quantum and classical evolution are identical in form, the classical evolution may already 
be started at time zero (r = —Mt c ), which is what we do in the numerical simulations 
(this is the analog of the "just the half" initial conditions in the case of the instantaneous 
quench @). 

An important point to make is that energy is not conserved during a mass quench of 
the type (|2.3j), because /i e fr depends explicitly on time. Differentiating the Hamiltonian 
corresponding to the action (2.1) and using the equations of motion one finds 

dH /> x ^k^. ( 2 .15) 



dt / dt 
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Hence during the linear quench there is a change in energy density 

Ap(t) = f dt> <^)(f) „ _^! 4 f dt'd 3 k\f k (t')\ 2 , (2.16) 

*Q j tQ JO 

where we made the Gaussian approximation in the last step (4 is the number of real 
components of the Higgs doublet). Evaluation of this analytic expression shows that the 
depletion in energy can be substantial 1 . 

In the case when the mass term comes from a coupling to an inflaton field ( |2.1| ), the 
energy extracted from the gauge-Higgs system is transferred into inflaton energy. Had 
we included the dynamics of the inflaton itself, this energy would not be lost but would 
come back through additional reheating as the inflaton oscillates and eventually decays. 
Including a dynamical inflaton does however complicate the system somewhat, since it 
implies unknown parameters such as A CT ^ and other masses and couplings of a potential 

VM). 

We shall use the approximation ( [2. 3D in the present work to bypass the complications of 
the additional degree of freedom, because it is easy to implement and solve analytically at 
early times, and because it can represent a source of mass quenching other than a coupling 
to an inflaton. We are therefore not bound to a specific model. 

When calculating the final asymmetry, we will need to know the final reheating tem- 
perature. We will then specialise to the case of a Higgs-inflaton coupling, assume m a > raw 
and that all the initial energy is equipartitioned between the SM particles at a temperature 
T rch < m w . 

3. No CP violation 

In order to understand the underlying gauge and Higgs dynamics, it is useful to study the 
evolution of various observables in absence of CP violation. The CP-violation can then, at 
least at early times, be thought of as a small perturbation on this background. 



The numerical implementation is identical to the one introduced in [14, gj. In short, the 
action (|2.l| ) is discretised on a lattice, and the classical equations of motion solved in real 
time. This is done for a CP-symmetric ensemble of initial conditions, reproducing the Higgs 
field correlators of the vacuum state at t = 0, when the Higgs potential is V = e + fjptfi'cf) 
in the Gaussian approximation. The initial gauge field is set by = and we impose 
the Gauss constraint on the gauge momenta for a given random Higgs field background. 
Observables are averaged over the ensemble. 

The observables of interest (in a periodic volume L 3 ) are the Higgs field averaged over 
the volume, 

~P=~^ J d 3 x4t(f>=-^ J <Fx~Tr&$, $ = (ir 2 0*,0), (3.1) 



1 The divergence of the integral over momenta is taken care of by renormalization, which may be ap- 
proximated by including only the unstable modes k < M y/r. The upper boundary of the integration over 
time has to be before nonlinearities become important, which is typically before the end of the quench at 
t = tg, except for a very fast quench. 
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the volume-averaged magnetic field, 

W=^JcPxTrF ij F ij , (3.2) 

the distribution of Higgs winding number, 

iV w = -^f d 3 xe ijk Ti E/t (Sjjcr) C/t [ft , jj = * (3.3) 

24tt J ^|lr$t$ 

the width of the Chern-Simons number distribution 



rt 

A cs (i) = ([iV cs (t) - ^(O)] 2 ), N cs (t) - N cs (0) = J Q dtJ d 3 x ^TtF^F^, (3.4) 

and its time derivative, the Chern-Simons number diffusion rate or susceptibility, 

t* dA(t) 



1 



dt 



(3.5) 



In equilibrium, T is the sphaleron rate. 

The dynamics of a tachyonic electroweak phase transition has been studied in some 
detail in terms of these and related observables in Jl3| , |23| (with an inflaton), in the limit of 
an instantaneous quench in [14, [24], |l|, and in terms of suitably defined particle numbers in 



J25[j . The above observables are all global quantities (integrals over space). Since = 
initially, also iV cs (0) = 2 . The Chern-S imons number and the winding number can be 
changed by an integer through a "large" gauge transformation, but N cs — N w is gauge 
invariant. Classical vacuum configurations have N cs — N w = 0, and are gauge equivalent 
to N cs = N w = 0. 

We now motivate the study of one other observable, the distribution of $<f> over the 



volume, in particular its magnitude near (f)'(j) = 0. The study in [24| concentrated on local 
observables, such as Chern-Simons and winding-number densities, in an attempt to clarify 
how CP violation causes an asymmetry in the final Chern-Simons number. It was observed 
that the transition produced initially many centers with high winding-number density, 
dubbed 'half-knots' since their winding number in small balls is roughly ±1/2. They can 
only disappear or be created when the Higgs length a/ 0t0 becomes zero in their center. 
This happens at early times because of the rapid growth of long distance modes, and then 
the number of half-knots rapidly diminishes, but when the volume- averaged Higgs length 
has grown substantially, (j) 2 = 0(v 2 ), such zeros have become rare. Subsequently the Higgs 
field 'overshoots' the minimum of its potential, rolls back and moves again towards zero, 
and when c/> 2 goes through its first minimum, new zeros in the Higgs length occur, enabling 
the creation of new half-knots. In this way second, third, . . . , generation half-knots were 
observed [^4|] near the minima of the oscillating (f> 2 . Since the gauge field is initially very 
small, CP violation is ineffective in influencing the creation or annihilation of the first 
generation half-knots. However, by the time of the first minimum of (j) 2 , the gauge field has 



2 We recall that iV cs = - / d 3 x e jki Tr Aj (F k i + i^AkAi) /16tt 2 
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grown substantially, CP violation can be effective and may change the balance and cause 
an asymmetry in {N cs ). Because of the importance of Higgs zeros, we shall in section 4.4 
present results for the distribution of the local Higgs field length, a simple histogram of 
(jr(x)(j>{x) for various quench times. 

We found in that the asymmetry depends strongly on the Higgs mass. The extreme 
cases seem to be mH = 2my/ ~ 161 GeV and = \/2mw ~ 114 GeV 3 . In addition to 
the quench time dependence, we are interested in how this mass dependence comes about. 
Below, we will show plots for these two Higgs masses in parallel, to demonstrate the 
differences throughout. Since the physical Higgs masses, at least in the Minimal Standard 
Model are constrained to be within 114 and 200 GeV S, both options are (marginally) 
allowed. 



3.1 Higgs field 




Figure 1: The normalized Higgs expectation value 2(<f> 2 )/v 2 vs time, for different quench times 
(full lines). Dotted lines show v 2 {t)/v 2 correspondingly colour coded. Left: m# = 2mw, right: 
run = \/2mw ■ 



As the Higgs mass parameter changes from positive to negative, the naive Higgs ex- 
pectation value goes from for t < ig/2, to v 2 (t) = —fil s (t)/X up until t = tQ, after which 
it stays at v 2 = /i 2 /A. 

Figure || shows the evolution of the normalized Higgs average value (2cj) 2 ) jv 2 (cf. ( |3.1| )) 
for various quench times. Dotted lines show v 2 (t)/v 2 for the different tQ. For small quench 
times, the symmetry breaking time is determined by the time it takes for the dynamics to 
perform the rolling off the local maximum of the potential at <fi = 0. For mHtQ < 36, v 2 (t) 
reaches v 2 before (2cj) 2 (t)} reaches v 2 (t). For larger quench times, (2(j) 2 (t)) catches up with 
v 2 (t), and oscillates around it before settling around v 2 . In the limit of infinite quench time 
one would expect it to follow v 2 (t) closely, except for finite-temperature corrections. 

We notice that the amplitudes of the first maximum and minimum are quench-time 
dependent, and also mass dependent. Although the qualitative behaviour is the same, 
tuh = \p2m-w leads to a lower first Higgs minimum. 

3 In 1+1 dimensions the mass dependence is complicated P2|, which could also be the case here. 
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The evolution of the Higgs field determines the energy loss in (the first half of) 
Eq. ( 2.16j ) due to the time-dependent effective mass. We integrate the actual numerical 
(4> 2 ) up to time t Q to find that for t Q = (0, 9, 13.5, 18, 36, 72) Eq. (|T6|) predicts 



|Ap/p initia i| ^ 0, 0.06, 0.10, 0.19, 0.66, 0.84. 
By directly calculating the energy, we find 

|Ap/p initial | ~ 0, 0.04, 0.08, 0.17, 0.67, 0.83. 



(3.6) 



(3.7) 



Note that for the very slowest quenches, more than half the energy is lost. For these 
slow quenches, at t = tQ the field has already started its oscillation, and so the Gaussian 
approximation (the second half of Eq. ( |2.16 )) does not apply. Integrating the actual field 
evolution reproduces the energy depletion. 

3.2 Magnetic field 



1.2- 2 
<B > 



— "Vcr 

- - "Vtf 9 




80 



m H t 



Figure 2: The evolution of the magnetic field B 2 for different tQ. ran = 2 m<w (left), tuh = V^mw 
(right). 

As the Higgs field goes through its transition, energy is transfered to the gauge fields, 
the occupation numbers of which grow exponentially pl| . The gauge fields acquire energy 
very fast at first, and then slowly towards what will be an equipartitioned and thermalised 
final state, figure |2[ This later transfer is quench-time dependent, faster quenches lead 
to a faster transfer of energy. To some extent, this decrease in energy transfer is due to 
the depletion of energy resulting from the time-dependent /i e ff) and the relative growth 
B 2 /B 2 is in fact similar for different quench times. The dependence on the Higgs mass 
is also interesting. For ran = V2mw the transfer seems to be driven by the Higgs field 
oscillations, with the gauge field also oscillating even for rather late times. In particular 
the second maximum is much larger than for run = 2mw- 

3.3 Chern-Simons diffusion 

In the absence of CP-violation, the ensemble average of the Chern-Simons number itself is 
zero, Chern-Simons number being CP-odd. In our setup, the ensemble is strictly CP-even, 
so (N cs ) is also strictly zero. 
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m H t 



Figure 3: (N^ s ) for different quench times. m# = 2niw (left), m# = \f2m\y (right). 



We can also calculate the evolution of the width of the N cs distribution, which will grow 
as a result of the preheating of the gauge fields, similar to the 5-field. But also because of 
fluctuations at finite temperature (and out of equilibrium, at finite energy density) there 
is a non-zero diffusion rate of Chern-Simons number, T(t) = dA cs /dt (cf. ( |3.4| )). The rate 
F(t) enters in an estimate of the baryon asymmetry in section ||. 

Figure || shows the evolution of A cs for various quench times. There is a rapid growth, 
which in some cases appears to be further driven by the Higgs field oscillations. Eventu- 
ally A cs settles, in accordance with the fact that at the final emerging temperatures the 
equilibrium sphaleron rate is negligible. This is one of the central features of Cold Elec- 
troweak Baryogenesis, the generated asymmetry does not get further diluted by sphaleron 
transitions. Similar to the magnetic field, at least for m# = \/2mw, the Chern-Simons 
number seems to be driven by the Higgs oscillations. However, whereas the growth of the 
magnetic field has a monotonic dependence on the quench rate, the diffusion rate appears 
to have a more complicated dependence. 



4. Adding CP-violation 



For small enough values of the CP-violation parameter <5 cp , it was seen in that the baryon 
asymmetry is linear. The value of <5 cp required to reproduce the observed asymmetry, is 
comfortably within this linear range. We shall keep 5 cp = 1 throughout, the upper end of 
the range studied in in order to maximise the generated numerical signal. 

Because the CP-violation is small, we have the option to treat it as a perturbation 
on the CP-even background described in the previous section. Although we will include 
the CP-violating term completely in the dynamics below, we will first apply the early time 
approximations introduced in [^6], |j| [l^, 14] for the case of finite-time quenches. 



4.1 Initial rise 

For very early times, only very long wavelength modes have large occupation numbers, 
both in the Higgs and gauge fields, since modes only gradually become unstable, starting 
with the zero mode at time £q/2. Treating the CP-violation as a perturbation to the 
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CP-even evolution, and making a homogeneous approximation, we find for the average 
Chern-Simons number |2^] 

V25 cp (Lm H f gg) (^2) 
^ cs; 64vr 4 (l+c) 2 v 2 /2" 1 j 

The values of 5 2 (t) and 2 (t) are taken from the simulations, section ||. The constant c 
is extracted from the growth of the gauge field B 2 {t) oc e 2ct . The linearisation assumes 
exponential growth of the Higgs and gauge fields. This description will have to break 
down at some fairly early time, when the back-reaction of gauge and Higgs non-linear self- 
interaction becomes important. Figure |] shows the linear approximation compared to the 
full simulation for short and long quench times. The agreement is good during the initial 
exponential growth, but breaks down after about 5 and 10 units of mnt after // 2 ff has gone 
negative at t = ig/2, respectively for £q = and 36. 



4.2 Thermodynamic treatment 




ny ny 

Figure 5: The average Chern-Simons number in the full simulations (full line) and from the 
thermodynamic treatment (dashed line), mji = 2m\y (left), \/2mw (right). 
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Beyond the linear approximation, we can apply methods from non-equilibrium ther- 
modynamics [26, ||, 13, ^] to estimate the asymmetry. 

One can interpret the CP-violating term as a chemical potential for Chern-Simons 
number 4 (cf. (OO)): 



(4.2) 



Using the CP-even evolution of the diffusion rate Eq. ([^) and the Higgs average Eq. ( |3~1 
the average Chern-Simons number can then be estimated through 



(N cs )(t) = 7 ±- / dt'T(t')^ ch (t'), 

J-eS 70 



(4.3) 



where T e g was interpreted in Q as the effective temperature of the tachyonic modes. We 
will not elaborate here on such an interpretation, but merely observe that T e fj turns out to 
decrease roughly linearly with tQ, and that m# = y2mw gives much larger values, figure 



Figure [5] compares the result of 
Eq. ( |Q| ) to the full simulation. T eS 
is chosen to fit the first maximum 
of the full simulation. The approxi- 
mation nicely reproduces the change 
of sign of the asymmetry produced 
by the back-reaction. At later times, 
the approximation again breaks down. 
We will see that this is precisely the 
time when the Higgs field acquires a 
net winding number the dynam- 
ics of which can apparently not be 
described by a simple chemical po- 
tential with constant T e g . The effec- 
tive temperatures as a function of tQ 
are shown in figure []. 

Notice in figure || that the sign of the asymmetry at later times rant ~ 40 has changed 
again to positive (the sign of <5 cp ) in the case of mass ratio \/2, which is not captured by the 
thermodynamic treatment. In principle the latter might do better, since the oscillations in 
Hch(t) an d r(i) are correlated. In any case, replacing the diffusion rate by its time average 

m 




Figure 6: The effective temperature in units of 
rriH as extracted from the thermodynamical treatment. 
Squares: ran = 2mw, circles: mu = V^rnw- 



tmax 



dt'T(t')^ ch (t') 



r / dt'» ch (t>) = ^4^, 

Zm w 



(4.4) 



gives a sign of the asymmetry that is definitely equal to that of 5 cp , which may be wrong 
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Figure 7: The evolution of Chern-Simons number in time, uih — 2mw (left), uih — V^Wff 
(right). 
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Figure 8: The evolution of winding number in time, m# = 2mw (left), m# = \/2mw (right). 



4.3 Full simulation 

In order to capture the full dependence on quench time and Higgs mass, we need to include 
the CP-violation completely in the dynamics. Figure ^ shows the average Chern-Simons 
number for various quench times. Figure || is the corresponding winding number. We notice 
that the mass dependence found in jlj is robust, and not a pathology of an instantaneous 
quench. For ran = 2mw, the fastest quenches mutq = 0,9 lead to an asymmetry of 
opposite sign to 5 cp . For slower quenches, the noise dominates and we can only conclude 
that the final asymmetriy is consistent with zero. In contrast, for run = y/2m\yi the 
asymmetry has the same sign as 5 cp , and it is maximal for intermediate quench times 
^HtQ = 18. Recall also the maximal boosting of A cs in section || was seen at these quench 
times. In both cases, for mntQ = 36 and larger the asymmetry appears to vanish. This is 
all compiled in figure ^ which shows the final asymmetry versus quench time. 

4.4 Higgs field zeros 

We have observed in earlier work [14, 24] that the final asymmetry in (N cs ) can already be 



'Notice that one treats <jy <j> as an space-independent chemical potential. 
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Figure 10: A close-up of the early evolution of (<fi 2 ) (black), (N cs ) (red) and (Nw) (blue). Full 
lines are mu = \^2mw, dashed m# = 2mw- The quench time mntQ = 9. The inset is a further 
amplification around the initial winding number bump. 



seen at earlier times in (N w ), which may be expected from the fact that the temperature 
after the transition is low enough that sphaleron transitions are suppressed and the robust- 
ness of winding number under relatively small changes in the fields. The asymmetry in 
(iV w ) is induced by the CP violating terms in the equations of motion, which are very small 
during the first stages of the instability, as monitored roughly by 2 and B 2 . Somewhat 
later the asymmetry becomes visible in the initial rise and bouncing back of (iV cs ), and a 
little later also in (iV w ). The rise in (N w ) is much smaller than in (N cs ), presumably since 
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iV w can only change when there are zeros in the Higgs field, which are exceptional by that 
time. Still somewhat later, around the time {(j) 2 ) has its first minimum, (N w ) has grown 
substantially, as new (second generation) zeros appear in the Higgs field. 

This is illustrated in figure pi], which shows the evolution of (iV cs ), (N w ) and (eft 2 ) at 
early times. At the time of the initial bump of (N cs ) the winding asymmetry (iV w ) is still 
not visible, but a growth of a 'bump' is discernible by the time (N cs ) has crossed zero 
and reached a negative maximum. When {(ft 2 ) reaches its first minimum the asymmetry in 
(N w ) grows much faster, which we interpret as being caused by the asymmetric creation 
of second generation winding centers, made possible by zeros in the Higgs lengths. 
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Figure 11: Histogram of (j> 2 (x) over the lattice for m# — 2m\y- Colours correspond to quench 
times mntq — (black), 9 (red), 18 (green) and 36 (blue). The four graphs correspond to the first 
four minima of the Higgs oscillation in each case. 

To illustrate the presence or absence of zeros in the Higgs field we show histograms 



of 4> {x) over the lattice for the first 4 minima of the Higgs oscillation, in figures 11 and 



|i2| . Figure |ll| is for ran = 2mw, and we see that although the average is away from zero 
(figure |l|), there is still a tail stretching to zero, at least for the first two minima. These will 
provide nucleation points for winding. In figure |l], we also saw that for the smaller Higgs 
mass ran = V^mw: the Higgs minima were somewhat lower. It is remarkable, however 



how different the distributions in the first minimum look (figure 12). The proximity of 
the distribution bulk to zero (and the fact that 4> 2 {x) > 0) results in points aggregating 
close to zero. Qualitatively, the density of zeros follows the same behaviour as a function 
of quench time as the final asymmetry of figure |9[ 

4.5 Kibble mechanism 

During a symmetry breaking transition, a net density of defects will form through the 
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Figure 12: Histogram of 4> 2 (x) over the lattice for m# = V2mw- Colours correspond to quench 
times mjjtg = (black), 9 (red), 18 (green) and 36 (blue). The four graphs correspond to the first 
four minima of the Higgs oscillation in each case. 
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Figure 13: The final distribution of winding number for various tQ. m# = 2mw- 



Kibble mechanism [27]. In the present case of 0(4) symmetry in 3+1 dimensions, these 
textures have integer winding number in the Higgs field, spread out over space. The density 
of gauged defects can be predicted in terms of the evolution of the correlation length of the 
system Numerical studies often consider thermal quenches with overdamped dynamics 
through the transition. In our case, we have an underdamped system with no explicit 



coupling to a thermal bath (higher-momentum modes do play the role of a bath) . Still, we 
here illustrate "fast" and "slow" quenches in terms of the number of defects generated. 

Figure 13 shows the distribution of final winding number over the ensemble. There 
is a qualitative difference between mutQ = 0, 9, 18 and mjjtQ = 36, 72. Also from this 
distribution we see that mutQ = 18 still belongs to the regime of 'fast' quenches, whereas 
miitQ = 36 is definitely in the 'slow' regime. 



5. Conclusions 

We have confirmed in a more realistic setting that including CP-violation in the SU(2)- 
Higgs equations leads to baryogenesis, when going through a tachyonic electroweak transi- 
tion. The dependence on quench time is significant. Although our scan of quench time is 
not fine enough to give a precise characterisation of 'fast' and 'slow' quenches, we can say 
that mutQ < 20 belongs to the first category, corresponding to 



df4s(t) 



2^dt 



> 0.07. (5.1) 

t=t Q /2 



This constrains a possible underlying Hybrid inflation model in terms of the Higgs-inflaton 
coupling A CT 0, since /^ ff (i) = fJ 2 - A^cr 2 . 

It is satisfying that the dramatic mass dependence seen in |l|] is not a result of the 
instantaneous quench. It is, however, surprising that for rriff = \^2mw the maximum 
asymmetry is not generated at the fastest quenches but for intermediate tQ. Both the 
mass dependence and the quench time dependence can be put down to a coincidence of 
phases and frequencies of the Higgs and gauge oscillations. This is similar to the much 
simpler model in 1+1 dimensions studied in pSfl . 

An important aspect of the transition is the occurrence of zeros in the Higgs field. 
CP-violation generates asymmetries in the Chern-Simons number density, which prompts 
the winding number density to move along as well. The full winding number change is 
however only realised once the second generation Higgs zeros appear. Although some 
settling of winding and Chern-Simons number can occur at later generations, we find that 
the asymmetry is established at the first Higgs minimum. 

The quench time and the mass (through the Higgs self-coupling) influence the number 
of zeros, and so determine the magnitude of the asymmetry by allowing more half- knots 
to flip their winding number. The detailed dynamics are very complicated. 

The expansion of the Universe is negligible at electroweak-scale temperatures, but 
reheating may proceed differently in the presence of all the other Standard Model fields 



and in particular an oscillating inflaton [23]. Still, because the final asymmetry is largely 
determined during the first couple of Higgs oscillations, we expect that the results presented 
here are reasonably close to the complete result. 

The maximal asymmetry occurred in our simulation for mu = V% m W an d rriHtQ = 18. 
We can estimate the photon density by distributing the initial energy density over the 
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Standard Model degrees of freedom 5 , which gives: 

— = (0.20 ± 0.04) x 10~ 3 5 cp . (5.2) 

n 7 

about three times 6 the result at zero quench time This means that we require 5 cp = 
3 x 10~ 6 or larger to reproduce the observed baryon asymmetry. 
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